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Hyperbolic formulations of the equations of motion are essential technique for proving the well- 
posedness of the Cauchy problem of a system, and are also helpful for implementing stable long 
time evolution in numerical applications. We, here, present three kinds of hyperbolic systems in 
the Ashtekar formulation of general relativity for Lorentzian vacuum spacetime. We exhibit several 
(I) weakly hyperbolic, (II) diagonalizable hyperbolic, and (111) symmetric hyperbolic systems, with 
each their eigenvalues. We demonstrate that Ashtekar's original equations form a weakly hyperbolic 
system. We discuss how gauge conditions and reality conditions are constrained during each step 
toward constructing a symmetric hyperbolic system. 
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I. INTRODUCTION 

Developing hyperbolic formulations of the Einstein equation is growing into an important research areas in general 
relativity Q. These formulations are used in the analytic proof of the existence, uniqueness and stability (well- 
posedness) of the solutions of the Einstein equation ||]. So far, several first order hyperbolic formulations have been 
proposed; some of them arc flux conservative some of them are symmetrizablc or symmetric hyperbolic systems 
P-im . The recent interest in hyperbolic formulations arises from their application to numerical relativity. One of the 
most useful features is the existence of characteristic speeds in hyperbolic systems. We expect more stable evolutions 
and expect imprements boundary conditions in their numerical simulation. Some numerical tests have been reported 
along this direction fl^lT^ . 

Ashtekar's formulation of general relativity has many advantages. By using his special pair of variables, the 
constraint equations which appear in the theory become low-order polynomials, and the theory has the correct form 
for gauge theoretical interpretation. These features suggest the possibility for developping a nonperturbative quantum 
description of gravity. Classical applications of the Ashtekar's formulation have also been discussed by several authors. 
For example, we [ p^ discussed the reality conditions for the metric and triad and proposed a new set of variables for 
Lorentzian dynamics. We ||l7[| also showed an example of passing a degenerate point in 3-space by locally relaxing 
the reality condition. Although there is always a problem of reality conditions in applying Ashtekar formulation to 
dynamics, we think that this new approach is quite attractive, and broadens our possibilities to attack dynamical 
issues. 

A symmetric hyperbolic formulation of Ashtekar's variables was first developped by Iriondo, Leguizamon and Reula 
(ILR) . They use anti-Hermiticity of the principal symbol for defining their symmetric system. Unfortunately, in 
their first short paper ]l^, they did not discuss the consistency of their system with the reality conditions, which are 
crucial in the study of the Lorentzian dynamics using the Ashtekar variables. We considered this point in [p9| , and 
found that there are strict reality constraints (alternatively they can be interpreted as gauge conditions). Note that 
we primarily use the Hermiticity of the characteristic matrix to define a symmetric hyperbolic system, which we think 
the more conventional notation. The difference between these definitions of symmetric hyperbolicity is commented in 
Appendix C. 

The dynamical equations in the Ashtekar formulation of general relativity are themselves quite close to providing 



a hyperbolic formulation. As we will show in §[V, the original set of equations of motion is a first-order (weakly) 
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hyperbolic system. One of the purposes of this paper is to develop several hyperbolic systems based on the Ashtekar 
formulation for Lorentzian vacuum spacetime, and discuss how gauge conditions and reality conditions are to be 
implemented. We categorize hyperbolic systems into three classes: (I) weakly hyperbolic (system has all real eigenval- 
ues), (II) diagonalizable hyperbolic (characteristic matrix is diagonalizable) , and (III) symmetric hyperbolic system. 
These three classes have the relation (III) e (II) e (I), and are defined in detail in §||. As far as we know, only a sym- 
metric hyperbolic systems provide a fully well-posed initial value formulation of partial differential equations systems. 
However, there are two reasons to consider the two other classes of hyperbolic systems, (I) and (II). First, as we found 
in our previous short paper [p^ , the symmetric hyperbolic system we obtained using Ashtekar's variables has strict 
restrictions on the gauge conditions, while the original Ashtekar equations constitute a weakly hyperbolic system. We 
are interested in these differences, and show how additional constraints appear during the steps toward constructing 
a symmetric hyperbolic system. Second, many numerical experiments show that there are several advantages if we 
apply a certain form of hyperbolic formulation. Therefore, we think that presenting these three hyperbolic systems 
is valuable to stimulate the studies in this field. To aid in possibly applying these systems in numerical applications, 
we present characteristic speeds of each system we construct. 

ILR, in their second paper [po| , expand their previous discussion concerning reality conditions during evolution. 
They demand that the metric is real- valued (metric reality condition) , and use the freedom of internal rotation during 
the time evolution to set up their soldering form so that it forms an anti-Hermitian principle symbol, which is 
their basis to characterize the system symmetric. However, we adopt the view that re-defining inner product of 
the fundamental variables introduces additional complications. In our procedure, we first fix the inner product to 
construct a symmetric hyperbolic system. As we will describe in our symmetric hyperbolic system then requires 
a reality condition on the triad (triad reality condition), and in order to be consistent with its secondary condition 
we need to impose further gauge conditions. The lack of these constraints in ILR, we believe, comes from their 
incomplete treatment of a new gauge freedom, so-called triad lapse A'^ (discussed in § HI ) , for dynamical evolutions in 
the Ashtekar formulation. In Appendix we show that ILR's proposal to use internal rotation to re-set triad reality 
does not work if we adopt our conventional definition of hyperbolicity. 

The layout of this paper is as follows: In we define the three kinds of hyperbolic systems which are considered in 
this paper. In §^j_we briefly review Ashtekar's formulation and the way of handling reality conditions. The following 
sections §[V and §0 are devoted to constructing hyperbolic systems. Summary and discussion are in j jVl[ Appendix 
^ supplements our proof of the uniqueness of our symmetric hyperbolic system. Appendices |^ and |^ are comments 
on ILR's treatment of the reality conditions. 



II. THREE DEFINITIONS OF HYPERBOLIC SYSTEMS 



We start by defining the hyperbolic systems which are used in this paper. 

Definition 1 We assume a certain set of (complex) variables Ua {a = l,---,n) forms a first-order (quasi-linear) 
partial differential equation system, 

dtUa,=f^^{u)dlUfJ+K^{u), (2.1) 

where J (th e ch aracteristic matrix) and K are functions of u hut do not include any derivatives of u. We say that 
the system (2.1) is: 

(I), weakly hyperbolic, if all the eigenvalues of the characteristic matrix are real pl| ]. 

(II). diagonalizable iiyperbolic, if the characteristic matrix is diagonalizable and has all real eigenvalues p^. 
(Ill), symmetric hyperbolic, if the characteristic matrix is a Hermitian matrix 

Here we state each definition more concretely. We treat J'^^q, as a n x n matrix when the /-index is fixed. The 
following properties of these matrices are for every basis of Z-index. 

We say A' is an eigenvalue of J''^q, when the characteristic equation, det ( J^'^q — \^5^a) — 0, is satisfied. The 
eigenvectors, p", are given by solving = A'p'". 

The weakly hyperbolic system, (I), is obtained when has real spectrum ioi every I, that is, when this characteristic 
equation can be divided by n real first-degree factors. For any single equation system, the Cauchy problem under 
weak hyperbolicity is not, in general, C°° well-posed, while it is solvable in the class of the real analytic functions and 
in some suitable Gevrey classes, provided that the coefficients of the principal part are sufficiently smooth. 

The diagonalizable hyperbolic system, |(II)| , is obtained when J is real diagonalizable, that is, when there exists 
complex regular matrix P' such that ((/") )"-y ,P'^s P^^ t3 is real diagonal matrix for every I. We can construct 
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characteristic curves if the system is in this class. This system is often used as a model in the studies of well- 
posedness in coupled linear hyperbolic system. (This is the same as strongly hyperbolic system as defined by some 
authors |^,^, but we use the word diagonalizable since there exist other definitions for strongly hyperbolic systems 



In order to define the symmetric hyperbolic system, (III), we need to declare an inner product {u\u) to judge 



whether J'^^q, is Hermitian. In other words, we are required to define the way of lowering the index a of u". We 
say J^^ a is Hermitian with respect to this index rule, when pa — J'' ap for every Z, where the overhead bar denotes 

complex conjugate. 

Any Hermitian matrix is real diagonalizable, so that |(III)| G |(II)| e |(I)|. There are other definitions of hyperbolicity; 
such as strictly hyperbolic or effectively hyperbolic, if all eigenvalues of the characteristic matrix are real and distinct 
(and non-zero for the latter). These defin ition s are stronger than |(II)| , but exhibit no inclusion relation with |(III)| . In 
this paper, however, we only consider |l)-(III) above. 



The symmetric system gives us the energy integral inequalities, which are the primary tools for studying well- 
posedness of the system. As was discussed by Geroch [^, most physical systems can be expressed as symmetric 
hyperbolic systems. 

III. ASHTEKAR FORMULATION 
A. Variables and Equations 



The key feature of Ashtekar's formulation of general relativity |15| is the introduction of a self-dual connection as 
one of the basic dynamical variables. Let us write the metric using the tetrad E'^, with i?^ satisfying the gauge 
condition i?" = 0. Define its inverse, E^, by g^^ ~ E^^E'lrju and i?j :— E^^gf^^rju, where we use /i, ^ = 0, • ■ • , 3 
and i,j = 1, • • • , 3 as spacetime indices, while I,J— (0), • ■ ■ , (3) and a,b = (1), • • ■ , (3) are 5*0(1, 3), SO{3) indices 
respectively. We raise and lower ^, j/, • • • by gf^" and g^i/ (the Lorentzian metric); I,J,--- by ij^'-' = diag(— 1, 1, 1, 1) 
and rjij; i,j, • • • by 7*-' and (the 3-metric); a,b, - ■ ■ by (5°^ and dab- We also use volume forms eabc- ^abc^""^^ — 3!. 
We define S0(3,C) self-dual and anti self-dual connections 

^A1■.^u';^T{^me\,u\^, (3.1) 

where is a spin connection 1-form (Ricci connection), Loj/ := E^^V ^E^ . Ashtekar's plan is to use only a self-dual 
part of the connection and to use its spatial part ^A°; as a dynamical variable. Hereafter, we simply denote 
as Al- 

The lapse function, N , and shift vector, A^*, both of which we treat as real-valued functions, are expressed as E^ — 
{1/N, —N^/N). This allows us to think of Eq as a normal vector field to E spanned by the condition t — x'~' =const., 
which plays the same role as that of Arnowitt-Deser-Misner (ADM) formulation. Ashtekar treated the set {E^, Af) 
as basic dynamical variables, where i?* is an inverse of the densitized triad defined by 

K ■■= (3.2) 

where e := deti?f is a density. This pair forms the canonical set. 

In the case of pure gravitational spacetime, the Hilbert action takes the form 

S = J d^xiidtAfWa + {i/2)NElElF^^e'^\ - e^AN - N^F^Ei + A-^V^K]. (3-3) 

where N er^N, F^^^ {dA"),,,, - (V2)e"bc(^'' A ^^)^^ = df.A'^ - d^A^ - ie'^bcA^.Al is the curvature 2-form, 

A is the cosmological constant, ViEl diEl - itab'^ A\El, and — detE* — (det-Ef)^ is defined to be detE* 
{l/6)e°''"'eijkElElE^, where eijk := CabcEfE^E^ and e^fe := er'^eijk [When {i,j,k) = (1,2,3), we have e^fe = e, 

e.jk: = 1, e'^'' = e-\ and = 1.]. 

Varying the action with respect to the non-dynamical variables N, N'^ and Aq yields the constraint equations, 

Ch ■■= (^/2)6''^ KElF^j - A det E « 0, (3.4) 
CAf. := -E^E^ « 0, (3.5) 
Caa := V,El « 0. (3.6) 
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The equations of motion for the dynamical variables {E^ and Af) are 

dtK = -iT^jie^'a NEiEl) + 2V,{N^^E'}) + i^tab^ El, (3.7) 
dtAf = -ie'^'^c NElF^j + N^F!^^ + + KNE^, (3.8) 

where VjXi' d^Xi' - ie^b'^^JXf , for X^^ + = 0. 

In order to construct metric variables from the variables {El,Ai,N,N'^), we first prepare tetrad Ej as E^ — 

(l/eTV, -NyeN) and E^ = (0,^^/e). Using them, we obtain metric gf'" such that 

g^'' := iJj'S}?]^'^. (3.9) 



B. Reality conditions 



The metric (3.9), in general, is not real- valued in the Ashtekar formulation. To ensure that the metric is real-valued, 
we need to impose real lapse and shift vectors together with two conditions (the metric reality condition); 



Im(^:^-'") = 0, 
Re(e"''^^^^f 2?fc^^)) = 0, 

where the latter comes from the secondary condition of reality of the metric Im{dt{El&°')} 
detE > (see 



J^). These metric reality conditions, (3.1C) and (|3.1l| ) 
evolution if the variables satisfy the conditions on the initial data |p8| , |16| . 

For later convenience, we also prepare stronger reality conditions. These conditions are 

i^K) = 

and Im(at£;^) = 0, 



(3.10) 
(3.11) 

[|8| , and we assume 
are automatically preserved during the 



(3.12) 
(3.13) 



and we call them the "primary triad reality condition" and t he " secondary triad reality condition " , res pe ctive ly. Using 
the equations of motion of E^, the ga uge c onstraint (3.4)-( ^^ , the metric reality conditions ( |3.10 ), (3.11) and the 
primary condition ( 3.1^ ) , we see that ( 3.13| ) is equivalent to [|16[| 



Re{A^) = di{N)E''' + {l/2e)E^NE^''djEi + N'Re{A'}), 



(3.14) 



or with un-dcnsitized variables, 



Re{A^) = di{N)E"' + N' Re(^^). 



(3.15) 



From this expression w e see that the secondary triad reality condition restricts the three components of "triad lapse" 
vector Aq. Therefore ( 3.14 ) is not a restriction on the dynamical variables {E^ and A";) but on the slicing, which 
we should impose on each hypersurface. Thus the secondary triad reality condition does not restrict the dynamical 
variables any further than the secondary metric condition does. 

Throughout this paper, we basically impose metric reality condition. We assume that initial da ta of {El , Af) fo r 
evolution are solved so as to satisfy all three constraint equations and metric reality condition ( 3.10| ) and (3.11). 
Practically, this is obtained, for example, by solving ADM constraints and by transforming a set of initial data to 
Ashtekar's notation. 



C. Characteristic matrix 



We shall see how the definitions of hyperbolic systems in §|| can be applied for Ashtekar's equations of motion (3.7) 
and ( |3.8| ). Since both dynamical variables, E^ and A'^, have 9 components each (spatial index: i = 1, 2, 3 and S0(3) 
index: a = (1), (2), (3)), the combined set of variables, — {E ljAf), has 18 components. Ashtekar's formulation 
itself is in the first-order (quasi- linear) form in the sense of (2.1), but is not in a symmetric hyperbolic form. 
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We start by writing the principal part of the Ashtekar's evolution equations as 



'K' 









\l hi 
L a 
f^lab 



El 
-45 



(3.16) 



where = means that we have extracted only the terms which appear in the principal part of the system. We name 
these components as A,B,C and D for later convenience. 
The characteristic equation becomes 



det 



bi 



O i'j 



= 0. 



(3.17) 



If ab'' and C'^-^ij vanish, then the characteristic matrix is diagonalizable if A and D are diagonalizable, since the 



spectrum of the characteristic matrix is composed of those of A and D. The eigenvectors for every ^-index, {p^]^,q^°'^ 
are given by 



A'a'^'j B'ab'' 



P"b 



la 



for every /. 



(3.18) 



The lowering rule for the a of u" follows those of the spacetime or internal indices. The corresponding inner product 
takes the form {u\u) :— Uau". According to this rule, we say the characteristic matrix is a Hermitian when 



— A''"'''^^ — A'^"-** 

Q j^labij ^Ibaji 



(3.19) 
(3.20) 
(3.21) 



IV. CONSTRUCTING HYPERBOLIC SYSTEMS WITH ORIGINAL EQUATIONS OF MOTION 



In this section, we consider which form of hyperbolicity applies to the original equat ions of motion, (3.7) and (^ 
under the metric reality condition (j jiV A ) or under the triad reality condition (§[VB). 



A. Under metric reality condition (system la and Ila) 

As the first approach, we take the equations of motion (3/7) and ( |3.8|) with metric reality conditions ( 3.1C ) and 
(3.11). The principal term of (3.7) and (|3.8D become 



dtAt = -le"^ NElF,'^ + N^Ff^ + + A NE^ 

^ -le'^K NElid^A'j - djA^) + N^idjA^ - O^A]) 



The principal terms in the notation of (3.16) become 



B 



dabij 



dabij 



0, 



We get the 18 eigenvalues of the characteristic matrix, all of which are independent of the choice of triad: 



(4.1) 

(4.2) 
(4.3) 
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(multiplicity=6), N' (4), N^±Ny^ (4 each), 



where we do not take the sum in 7" (and we maintain this notation hereafter for eigenvalues and related discussions) 



Therefore we can say that this system is weakly hyperbolic, of type [I] 



We note that this system is not type |(II)| in general, because this is not diagonalizable, for example, when iV' 
We classify this system as type (I), and call this system la, hereafter. 



0. 



The necessary and sufficient conditions to make this system diagonalizable, type (II) , are that the gauge conditions 



iVV nor ± Ny^r 



and 7" > 



(4.4) 



where the last one is the positive definiteness of 7". This can be proved as follows. Suppose that (4.4) is satisfied. 
Then 0, A^', A^' ± Vt"^ ^^"^ four distinct eigenvalues and we see rank(j') = 12, rank(j' — iV'/) = 14, rank(j' — 
(iV' ± N\/^)I) — 14. Therefore the characteristic matrix is diagonalizable. Conversely suppose that iV' = or 
TV' = iiVv^T", then we see the characteristic matrix is not diagonalizable in each case. 

The components of the characteristic matrix are the same as system la, so all eigenvalues are equivalent wit h sys tem 
la. We can also show that this syst em is not Hermitian hyperbolic. Therefore we classify the system [la + (|4.4| )] to 
real diagonalizable hyperbolic, type (II), and call this set as system Ila. However, we will show in the next section 
that real diagonalizable hyperbolic system can also be constructed with less strict gauge conditions by modifying 
right-hand-side of equations of motion {system lib). 



B. Under triad reality condition (system lb) 



Next, we consider systems of the original equati ons o f mo tion, (3.7) and (3.8), with the triad reality condition. 
Since this reality condition requires the additional ( 3.14 ) or ( p.lq ) as the secondary condition (that is, to preserve 
the reality of triad during time evolution), in order to be consistent with this requirement and to avoid the system 
becoming second order in fundamental variables, we need to set diN = 0. This fixes the real part of the triad lapse 
gauge as Re(ylg) — Re(^fiV'). We naturally define its imaginary part as Im(^g) = lm{Af N"^) . Thus the triad lapse 
is fixed as Aq — A^^N"^. This gauge restriction does not affect principal part of the evolution equation for E^, but 
requires us to add the term 

P.^S - d^A^o = ^^{A]N^) - N^d^A';) = N^S^SiidiA';) 



to the right-hand-side of the equation of A^. That is, we need to add N^S^Sl to I?'°bi^ in (4.3), 



D 



labij 



ibc 



The other components of the characteristic matrix remain the same [(4.1) and (4.2)]. We find that the set of eigenvalues 
of this system is 



(multiphcity =3), iV' (7), N' ± NVj" (4 each). 



Therefore the system is again, type (I). This system is not real diagonalizable because is not. So we classify this 
system as type [l]] and call this set as system lb. We note that this system is not real diagonalizable for any choice 
of gauge. Therefore we cannot construct a system of type (II) using the same technique of constructing system Ila. 
However, as we will show in the next section, the system becomes diagonalizable (and symmetric) hyperbolic under 
the triad reality condition if we modify the equations of motion. 



V. CONSTRUCTING A SYMMETRIC HYPERBOLIC SYSTEM 



From the analysis of the previous section, we found that the original set of equations of motion in th e As htekar 
formulation constitute a weakly hyperbolic system, type |[lj, or a diagonalizable hyperbolic system, type |(II)| , under 
appropriate gauge conditions, but we also found that we could not obtain a symmetric hyperbolic system, type |(III)| . 



In this section, we show that type (HI) is obtained if we modify the equations of motion. We begin by describing our 
approach without considering reality conditions, but we will soon show that the triad reality condition is required for 
making the characteristic matrix Hermitian. 
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We first prepare the constraints (p.4|)-(|3.6|) as 



CMk - -F^M = -{duA'; - d,Al)Ei = [SiEi + diEi]{diA% 

CGa^V.El^diEl 



(5.1) 
(5.2) 
(5.3) 



We apply the same technique as used by ILR to modify the equation of motion of and A"; ; by adding the constraints 
which weakly produce Ch ~ 0, Cmu ~ 0, and Cca ~ 0. (Indeed, this technique has also been used for constructing 
symmetric hyperbolic systems for the original Einstein equations |^,^.) We also assume the triad lapse Aq is 



- T^'^j diEl + S'\,^ diA' 



(5.4) 



where T and S are parameters which do not include derivatives of the fundamental variables. This assumption is 
general for our purpose of studying the principal part of the system. 

One n atural way to construct a symmetric hyperbolic s yste m is to keep B = C = and modify the A and D terms 
in ( 3.1(j ), so that we modify (3.7) using Co, and modify (3.8) using Ch and Cm- That is, we add the following terms 
to the equations of motion: 



modifying term for dtK = P\bC'a = P'a'djEf, ^ {P\'' 5]){diEi), 
modifying term for dtAI = ViAo + QfCn + Rt'" Cmj 

= T'''\j diEi + S'\^ diA'^ - tQ^eb'"' EIE'M^]) + + SiEi]diA 



(5.5) 



[S'\^ - iQW" EiE\ - i?> El + Rr Ei\{diA]) + T^'Kj diEl (5.6) 



A2 



we show that the modifications to the 



where P, Q, and R are parameters and will be fixed later. In Appendix 
off-diagonal blocks B and C, i.e. modifying (3.7) using Ch and Cm and modify ( |3.8| ) using Cg, will not affect the final 
conclusion at all. Note that we truncated Aq in (5_^), while it remains in (|^), since only the derivative of Aq effects 
the principal part of the system. The terms in (3.16) become 



AlaUj ^ ^i^bcajy^ljj^i^ _ ^^cba ^ j^l^^rj ^ J^l^Hgab _ J^i^U gab _^ piab^lj ^ 

j^labij _ Q 
^labij rplabij 

jjlabzj ^ i^abcj^^^U _ i^abcj^ji^l^^n ^ J^^lgab^ji _ J^jgab^H _^ glabij 

-iQ'''e'""^EiEli - R'^'^E^'' + W'&K 
The condition (|3.2l|) immediately shows T'"''*-'' = 0. The condition (|3.19| ) is written as 
^ -ie^^^N^^^El + ie'''"=N-f''Ei - 2e'''"=Nj'^ Im(£;') 



(5.7) 

(5.8) 
(5.9) 

(5.10) 



(5.11) 



By contracting f'afc'^, we get Re(eQ6c f'"''*''' 7;,; — 2€abc ]^°-^^^ j^j) = 20 7yim(_E^). This suggests that we should impose 

Im(£'') — 0, in order to get f'afc*"' = 0. This means that the triad reality condition is required for making the 
characteristic matrix Hermitian. 



A. Under triad reality condition (system Ilia) 
In this subsection, we assume the triad reality condition hereafter. In order to be consistent with the secondary 



triad reality condition (3.15) during time evolution, and in order to avoid the system becoming second order, we need 



to specify the lapse function as dtN = 0. This lapse condition reduces to 



Re{A^)^N'Re{AT), (5.12) 
a, Re(^S) = N'd^ Re(^p. (5.13) 
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By comparing these with the real and imaginary components of ( |5.4| ), i.e., 

di Re(^g) ^ ReiS'^u') di Re{A)) - Im(5'%,^) di lui{A)), (5.14) 

ajm(^g) = lYn{S'\^)diRe{A])+Re{S'\^)dilu,{A]), (5.15) 

we obtain 

ReiS^'^M^) = N^Sl 6^ and Im(S''%,^) = 0. 

Thus S is determined as 

S'\^=N^lS^. (5.16) 

This value of 5, and T = 0, determine the form of the triad lapse as 

Aq — A°N' + non-dynamical terms. (5-17) 

ILR do not discuss consistency of their system with the reality condition (especially with the secondary reality 
condition). However, since ILR assume Aq = AfN"^, we think that ILR also need to impose a similar restricted lapse 
condition in order to preserve reality of their system. 

By decomposing |, that is ( 3.1£ ), into its real and imaginary parts, we get 

= -N'-f^^S""^ + N^-^^'S^" + 7'^' Re(P)''''^ - 7'* Re(P)^'''", 

= -e^'^'^N-i^^El - e'^'^Nj^'El + 7'^' Im(P)*'''' + 7'* Im(P)^''"'. 

By multiplying 7;,; to these and taking symmetric and anti-symmetric components on the indices ab, we have 
= 2N^^'"'^ + Re(P)^'("'') - 3Re(F)^('"^) = 2N^'"' - 2Re(P)^'("''\ 
= 2Ar^'(5['"^l + Re(P)^'[''^l - 3Re(P)^'['"'l = 4Re(P)^'[''''l , 
= Im(P)^("'') + 3Im(P)^'(''") = 4Im(P)^(''^\ 

= -2e'"'^NEi + Im(P)^["''l + 3Im(P)j[''"l = -2e'"'^NEi - 2Im(P)j[''''l. 

These imply 

Our task is finished when we specify the parameters Q and R. By substituting (5.16) into ( 5.10| ), the condition 
(3.20) becomes 

= ie'''"'NEl-i^' + ie'""'NEl-/^^ - iQ"' t^""^ E^E^ - iQ''^ e""'^ El&a - IT^^E^^ + R'^'^E^^ 

+W^''E"' - R^'''^". (5.19) 
We found that a combination of the choice 

g" = e-^NE"", and P^'" = ie~^ Ne'"='^EiiEi, (5.20) 



satisfies the condition ( 5.19 ). We show in Appendix A_l that this pair of Q and R satisfies ( 5.1£ ) and that this choice 
is unique. 

The final equations of motion are 

B'"'"^ = C'"''*^' = 0, (5.22) 
D'"'"^ = iN{e'''"'El-f'' - e''''=P^7^'* 

_g-2|;»ag6cd^j^^ - e-'^e'""^EiiElEi'' + e''^ e"""^ E'^Ei&^) + N^S'^'^Y'- (5.23) 

To summarize, we obtain a symmetric hyperbolic system, type KlII)| by modifying the equations of motion, restricting 
the gauge to: Aq = A'^N'^, diN = 0, and assuming the triad reality condition. We name this set system Ilia. The 
eigenvalues of this system are 



(multiplicity = 6), ± Vt"^^ (5 each) and ± a/t"^ (1 each). (5.24) 

These speeds are again independent of the choice of (real) triad. 
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B. Under metric reality condition (system lib) 



Using this technique, we can also construct another example of diagonalizable hyperbolic system. Since the param- 
eters iS* and T specify triad lapse, a gauge variable for time evolutions, it is possible to change our interpretation that 
we take the evolution of the system within the metric reality condition. Of course, the characteristic matrix is no 
longer Hermitian. From the fact that we do n ot use the triad reality condition in the process of modifying the charac- 
teristic matrix using parameters ( 5.18 ), ( |3.20 ) nor in the process of deriving the eigenvalues, this system has the same 
components in its characteristic matrix and has the same eigenvalues. The process of examining diagonalizability is 
independent of the reality conditions. Therefore this system is classified as a diagonalizable hyperbolic system, type 



(11) 



To summarize, we gain another diagonalizable hyperbolic system by modifying the equations of motion using terms 
from con straint equations, with characteristic matrix (5.21)-( 5.23| ) under metric reality condition. The eigenvalues 
are ( ^.24 ), and this system is restricted only by a condition on triad lapse, A'^ = AfN^, and not on lapse and shift 
vector like system Ila. We call this system system lib. 



VI. DISCUSSION 



We have constructed several hyperbolic systems based on the Ashtekar formulation of general relativity, together 
with discussions of the required gauge conditions and reality conditions. We summarize their features in Table |[ 



system 


Eqs of 


reality 


gauge conditions 


first 


all real 


diagonal- 


sym. 




motion 


condition 


required 


order 


eigenvals 


izable 


matrix 


la 


original 


metric 




yes 


yes 


no 


no 


lb 


original 


triad 


A?i = AtN\ ^^N = 


yes 


yes 


no 


no 


Ila 


original 


metric 


N' / 0, ±N^/ji^ (7" / 0) 


yes 


yes 


yes 


no 


lib 


modified 


metric 


A^ = AtN' 


yes 


yes 


yes 


no 


Ilia 


modified 


triad 


Ao=A'^N\ ^^N^0 


yes 


yes 


yes 


yes 



TABLE I. List of obtained hyperbolic systems. The system /, //and ///denote weakly hyperbolic, diagonalizable hyperbolic 
and symmetric hyperbolic systems, respectively. 



The original dynamical equations in the Ashtekar formulation are classified as a weakly hyperbolic system. If we 
further assume a set of gauge conditions or reality conditions or both, then the system can be either a diagonalizable 
or a symmetric hyperbolic system. We think such a restriction process helps in understanding the structure of this 
dynamical system, and also that of the original Einstein equations. From the point of view of numerical applications, 
weakly and diagonalizable hyperbolic systems are still good candidates to describe the spacetime dynamics since they 
have much more gauge freedom than the obtained symmetric hyperbolic system. 

The symmetric hyperbolic system we obtained, is constructed by modifying the right-hand-side of the dynamical 
equations using appropriate combinations of the constraint equations. This is a modification of somewhat popular 
technique used also by Iriondo, Leguizamon and Reula. We exhibited the process of determining coefficients, showing 
how uniquely they are determined (cf Appendix |A|) . In result, this symmetric hyperbolic formulation requires a triad 
reality condition, which we suspect that Iriondo et al implicitly assumed in their system. As we demonstrated in 
j |v| , in order to keep the system first order, and to be consistent with the secondary triad reality condition, the lapse 
function is strongly restricted in form; it must be constant. The shift vectors and triad lapse Aq should have the 
relation (5.17). This can be interpreted as the shift being fre e an d the triad lapse determined. This gauge restriction 
sounds tight, but this arises from our general assumption of (5.4). ILR propose to use the internal rotation to reduce 
this reality constraint, however this proposal does not work in our notation (see Appendices ^ and 

There might be a possibility to improve the situation by renormalizing the shift and triad lapse terms into the 
left-hand-side of the equations of motion like the case of general relativity Or this might be because our system 
is constituted by Ashtekar's original variables. We are now trying to relax this gauge restriction and/or to simplify 
the characteristic speeds by other gauge choices and also by introducing new dynamical variables. This effort will be 
reported elsewhere. 



We thank John Baker for his useful comments on the initial draft. We thank Abhay Ashtekar for his comments 
on our draft. We also thank Matt Visscr for careful reading the manuscript. A part of this work was done when HS 
was at Dept. of Physics, Washington University, St. Louis, Missouri. HS was partially supported by NSF PHYS 
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96-00049, 96-00507, and NASA NCCS 5-153 when he was at WashlJ. HS was supported by the Japan Society for the 
Promotfon of Science. 



APPENDIX A: DETAIL PROCESSES OF DERIVING THE SYMMETRIC HYPERBOLIC SYSTEM IIIA 

In this Appendix, we show several detail calculations for obtaining the symmetric hyperbolic system Ilia. 



1. Determining Q and R 



We show here that the choice Q and R of (5.2C) satisfies (5.19). That is, the final D in ( p. 23 ) satisfies (3.2C), and 
that this choice Q and R is unique. 

First we show that D in ( ^.23| ) satisfies Hermiticity, ( ^3.20 ). From the direct calculation, we get 



e-h'""^E"'EiE'a 



'e^'^E^'ElE'a - e-^^^^K^ElEl) 



+iN{e'""'Ei-f^^ - e'""'Eif^ - e-''e'""'E^''ElE'a + e-'^""" E"-" E'^E'^ - e^^e'^^S'^SJi^J) 



— 2 ^bcd ^ia ] 



abc Ij 



-'-'^e'""^&^ElE\ - e-''e'"'''E'''ElE]i) ^-.iN] 



Llahij 



(This -j-'^^y definition is used only within this subsection Al,) Hermiticity, Y'^^'^-' = 0, can be shown from the fact 

2 Yiab)tj ^ -e-^€'""^&^ElEli - e~^e'""^E"'EiE'^ - e'^ e'""^ E^E'^ - e'"^ e'""' El&^ = 0, 
and its anti-symmetric part 0, which is derived from 

±labij abc Tpj li abc jpi Ij —2 acd jjilb j^i -rpj —2 bed jjila jjij j^i 

= 2Ei-f" - 2EIj'^ - e-^&^ElEi + e-^&^ElE'i - e-^E'^^E^El + e-^E'^'EiEl 



= 2EiY' 



Next we show that the choice Q and R of (5.20 ) is u nique in order to satisfy ( 5.19] ). Suppose we have two pairs of 
(Q, i?), say {Q i,Ri) and ((52, -R2), as solutions of ( 5.19| ). Then the pair {Qi — Q2, Ri — R2) should satisfy a truncated 
part of (pT9|), 



X 



labij 



; Q''^e'""^Ei&^ - i Q^^e'""^El&^ - R'^^'P'' + R'^''&'' + W^'^E"' - W = 0. 



Now we show that the equation f-"-^"^! — has only the trivial solution Q = R = 0. By preparing 

^abtj ^ -iQ^'e^^-^ElEd, - R'^-'-iuE^^ + R'^^E^, 



we get R^^°'^ij — 0. By substituting this into (A2), we can express i? by Q as 



-^ — 2/~\ak ^bcd 



e'-^EiEkdEl^iQ''''e^h 



(Al) 
(A2) 

(A3) 



Therefore (Al) becomes 



f"'"^ = -i Q'''e'""^Ei&a - i Q^^e'""^El&d ~ * Q^^^'^k E^^ + i Q"^?-^ k - i Q^'"e^\ E'" 



From this equation, we get the following contracted relations: 

^-2 ^labij ^ _i n^-i^bjl _ ^ ^bj ^ail _ ^ ^ak^il ^ pjb _^ ^ ^ak Jj ^ plb 

g-2 ^labtj ^ _ ■ qa^jl 4. 2i Q^^ea^'k ^^''^ " 2i Q^'^e^' 
^ — 2 xlabij TP ^_ . o^" /"iQ. ^b i o^" r^bc o^' r^cb i aa /Sbc 



i Q'^'e^'k E^" + i Q^'^e'^fc E'" - i Q'"'e^'k E'" + i Q^^t^\ E^" 



t'""*^ E^aeijc = 2i Q\S', + 2i Q"' ~ 2i Q'" + 4i Q" 



(A4) 



— 2 xlabij 



r"'' E,aeijc5t = &iQ\+AiQ\ - 10zRc(Q\) - 2Im(Q\). 



where Q°-^ := Q°-^E^ and e^^^ e^^^E^. From the last one, we get = 0. By substituting this into (A4), we get 
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-.be 



eb 



.[be] 



(A5) 



The symmetric part of (A5) indicates Q^^'^-' — 0, and 

g-2 ^labrj ^^^^^^ ^ 2i Q''=eb,' - 3i Q''^e6c' = -Re(Q^^ete') + 5i Im(Q^^efc,*) 
gives us Q^'^ebe^ — Q^''^' = 0. Therefore Q^^ — Q"' = is determined uniquely. From ([A^), we also get = 0. 



2. Modifications to ofF-diagonal blocks 



On the starting point of the modifications to the equations of motions (5.5) and ( |5.6| ), we assumed that off-diagonal 
terms keep vanishing. In this subsection, we show that the modifications to the off-diagonal blocks B and C in the 
matrix notation of ( p. 16 ), i.e. modifying (3.7) using Ch and Cm and modify (3.8) using Cq, does not affect the final 
conclusion at all. 



Suppose we have a sy mmetric hyperbolic system ( 5.21 )-( 5.23 ), and suppose we additionally modify the equations 
of motion (3.7) and ( |3.8| ) as 



modifying term for dtK = G\Ch + K'Cmj 



G\{~ie,^'' EiE\){diA]) + Hf{-5iEi + 5iEi){diA]) 



modifying term for dtA'', = I'^^Cab = (1"^ S]){diE^ 



bh 



where G^, H^J and I i are parameters to be determined. In the matrix notation, these can be written as 

B\b'' = -iG^a^,^' El&,-Hl'El+ Hl^El 



C 



lab 



jab xl 



The Hcrmitian condition (3.21) becomes 



= -iG],eb'^ M&d - K'El + WJEl - ha' 7'^ =: fab'' 



(A6) 
(A7) 



(A8) 
(A9) 



(AlO) 



(We use this Yab^' definition only inside of this subsection A 2.) 

If there exists a non-trivial combination of G^, H^^ and I'^^'i which satisfy this relation, then it will constitute 
alternative symmetric hyperbolic system. However, we see only the trivial solution is allowed for (AlO) as follows. From 
the relations of t'"'^*^7»j + t'°'"''7H = -iP"'', we obtain = 0. With this = 0, we obtain fab'' E] = -'2e^H'^\ 
which determine H^^ = 0. Similarly, from = and iJ*^ ~ 0, we get Ya''' ^jikE^ ~ ^Sie^G^, which determine 
Gi =0. 



APPENDIX B: INTERNAL ROTATION AND ASHTEKAR EQUATIONS 

In this Appendix, we consider the effect of a 50(3) rotation on the triad, which corresponds to a SU(2) rotation 
on the soldering form. The equations that we derive here will be applied in the discussion in Appendix O. 



1. Primary and secondary conditions of internal rotation 



The 50(3) internal transformation only affects inner space, and not the space-time quantities. Let us write U for 
such a rotation. U should satisfy the condition 



ab 



(Bl) 



This comes from the transformation of 6°-^ to S*°-^ := If^c U^d.S'^'^i which should satisfy 5*"'' = 5°''^. The determinant 
det U must be ±1, and we cho ose d et C/ = 1 for later convenience. The transformation 6°'b ~* S*°'b is naturally defined 
by S*''b ■■= WcUb'^S'^d- From we get the fundamental relations: S*''b = S^b, S*ab = Sab, and e*"''" = e"''". 
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Now we define the transformation of tfie triad El and of the inverse triad E? as 



E*f 



Ua'El 

U\El 



(B2) 
(B3) 



The S-metric, 7'^ , is preserved under this transformation, since 7'-' = E^^E^"- = E*\E*^°- . We note that this secondary 
condition, 9t7*-' — dt{El^E^°') ~ dt{E*l^E*^"-), wiU not give us further conditions. This is equivalent with the time 
derivative of (Bl). 



2. Internal rotation of Ashtekar variables 



Using det U = 1, the transformation of the densitized triad becomes 

EZ^Ua'El 

and straightforward calculation shows 

A*'} = U\A':--'-e'''',Ub'id,U^i), 



(B4) 



(B5) 



where we also note that a;*°° = U°-b and 0;*^'= = C/°e(e*'&c w,''^) - e''bc{dzU^'^)U''d- We remark that the second term 
in ( |B5|) arises because Af includes the spatial derivative of the triad. The relations of triad lapse and curvature 2-form 
become 



A*^o^U\A'o~-e'''cUb'idtU'=a), 



and constraints (3.4)-(3.6) are transformed into 



C* Mi = Cah, 
C*Ga — Ua^ Ccb- 



(B6) 
(B7) 



(B8) 
(B9) 
(BIO) 



The Hilbert action ( pT^ ) will be preserved (S* — S) under U, which is demonstrated by the "cancellation relation" 

idtA*t)E*l + C*Ga = idtAt)K + ^ Cca. (BU) 



Therefore the equations of motion for i?*^ and A*f are equivalent with the original ones, (3/7) and (3^), putting a * 

on all terms. 

The secondary metric reality condition ( |3.1l| ), W'^ := Re{e'''"'E*'^E*['V*kE*c), retains its form. 



while the secondary triad reality condition (3.15), := — Rc(^q) + di{N)E^"- + A^*Re(^"), is transformed as 



Y* 



Re{U''b)Y'' - ^^^{N)Re{U''b)lm{El) +lin{U\)[lin{Al) - ^^iN)lm{E''') - N' Im(A^)] 



1 



+ -A^'e%,Im(t/'"^)(9am(C/^d)). 



(B12) 



This equation has many unexpected terms, even if we assume the triad reality, Im(i?* ) = 0, before the transformation. 

To summarize, under triad transformations, Af, Aq, and are not transformed covariantly, while the other 
variables are transformed covariantly. 
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3. Make triad real using internal rotation 



Suppose all the variables satisfy the metric reality conditions, that is, satisfies Ini(£'* E'^") = 0. Can we obtain 
the triad which satisfies the triad reality condition, lm{E*l^) = 0, by an internal rotation? 
The answer is affirmative. However, such a rotation U must satisfy 

= Im(^*^) = lm{Ua^ El) = Re{Ua^) lm{El) + lm{Ua'') Re{El), (B13) 

and its secondary condition 

= lm{dtE"a) = lm[{dtUj)El + Ua' (dtEl)]. (B14) 



The application of this technique will be discussed in §|C2|. Before ending this section, we remark two points. First, 
Af is not transformed covariantly by this rotation U. Second, when we consider the evolution of E*l^, the evolution 



should be consistent with the secondary triad reality condition (3.14) 



APPENDIX C: CONSIDERATION OF ILR'S TREATMENT OF REALITY CONDITIONS 

The symmetric hyperbolic system (system Ilia) that we obtained in is strictly restricted by the triad reality 
condition. ILR (in their second paper [^ ) propose to use an internal rotation to de-constrain this situation. Here 
we comment on this possibility. 

1. Difference of definition of symmetric hyperbolic system 

First of all, we should point out again that there is a fundamental difference in the definition used to characterize 
the system as symmetric. As we discussed in §||, we define symmetry using the fact that the characteristic matrix is 
Hermitian, while ILR ||lj,|2^ define it when the principal symbol of the system iB'j°'ka {iJ^^ aki in our notation) is 
anti-Hermitian. 

We suspect that these two definitions arc equivalent when the vector ka {h in our notation) is arbitrary real. 
Actually, ILR have advanced a suggestion that our definition and their 'modern' version are equivalent. The judgement 
which is conventional or not, however, we would like to leave to the reader. Concerning our definition of symmetric 
hyperbolicity, we think that the readers can quite easily compare our system with other proposed symmetric hyperbolic 
systems in general relativity: all eigenvalues (in the system we presented) are all real- valued, while ILR's a re all pure 



imaginary. (Even if the distinction of real and pure imaginary is ignored, the eigenvalues calculated by us (5.24) and 
by ILR are different.) 

We note that, in addition, this fundamental difference will lead to different conclusions regarding the treatment of 
the reality condition (see the proceeding discussion). 

2. Can we obtain a symmetric hyperbolic system by internal rotation? 

What ILR proposed is the following: Suppose the system satisfies the reality condition on the metric, but not on the 
triad. By using the freedom of making an internal rotation, we can transform the soldering form to satisfy the triad 
reality condition, in such a way it forms symmetric hyperbolic system. (In their terminologies, they seek a "rotated" 
scalar product that is to find a more general symmetrizer.) Therefore we can remove the additional constraints of the 
triad reality. 

This procedure, however, includes changing inner product of dynamical variables, which might cause the topology 
of well-posedness of the initial value formulation to change. Here, we examine whether such a re-definition of the 
inner product is acceptable in our definition of symmetric hyperbolicity. 

Suppose we have a system which satisfies the constraints, and the metric reality condition, but not the triad reality 



conditions. As we commented in §[II, metric reality will be preserved automatically by the dynamical equations 
( |3.16D and (|5.21|)-(|5.23|). Now we apply a SO{3) rotation E^ E*\ := Uj' El to the system. We summarized 
the transformations of Ashtekar's variables and equations by U in Appendix |b[ In the new variables {E*\, A*1), 
transformed via J7, the equations of motions are written covariantly. 

As discussed in § p3 3| , it is possible to constru ct the real triad by using U . However , we a lways should verify the triad 
reality condition, both its primary condition ( |3.12 ), and its secondary condition ( 3.13| ). The latter is expressed as 
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(B12) or (B14). If we interpret this secondary condition as a restriction on the gauge variables, lapse N, shift N"^, and 
triad lapse ^q, then we only need to solve the primary condition in order to obtain triad reality on 3-hypersurface. 
This is indeed solvable. For example, ILR explain a way to get a real triad using orthonormality of the basis in their 
Appendix A in [^ . 

Next, let us see whether a symmetric hyperbolic system is obtained by the new pair of variables {E*l^, A*i). We 
define the equations of motion similarly as 



dt 







A*1 





A-^l hi 
^ a j 



d. 



terms with no diE*^ nor diA*f. 



(CI) 

(C2) 
(C3) 

(C4) 

These equations are related to (5.21)-( 5.23| ). We note that, in the modification here, we added the terms 
{N^S^'' + iNe^'"^E*l)CQi, coming from the terms of the gauge constraint. This corresponds to the relation A*'-"-''^^ — 

VcU^dA^^'^^i . 



By applying the same modifications as those in we get 

j^^labij ^^labij q 

jj*labij _ lJ\[(^^°-bcj^*j^^li _ ^abcj^*l^^ji 
-2 Tp^ia J)cd 



-e ■'E*'''e'""'E*lE ^ 



e ~e—E*i,E*iE*^''- 



e-^e'""^E"aE*iE*^'') + N'S'^''^^. 



Equations (|C2|)-(C4) forms a Hermitian matrix in the principal part of (CI), but it contradicts the consistent 
evolution with triad reality. That is, for example, the left-hand-side of dynamical equation dtE*\ = ■ • ■ [upper half 
of (|Cl|)] is real-valued since we impose Im(£^*^) = 0, while in the right-hand-side includes complex value in the non- 
principal part. To explain this in another words, the system (C1)-(C4) will not preserve the triad reality. Therefore 
we again need to control gauge variables through the secondary triad reality condition, and this discussion again 
returns the same gauge restrictions with those in § [v| . 

We also point out that the inner product of the fundamental variables in our notation does not form Hermitian like 
in the case of ILR. The inner product before the rotation U can be written 



which is common to ours and ILR's, while after the rotation the inner product becomes 

{iE*l,A*mE*lA*'})) = U,''U^'KEl + Uc''U^''l''A'^A'^ 



(C5) 



1 



1. 



(C6) 



which is not Hermitian, and can not be used as the inner product of the original variable (_E*,^") as in the ILR's 
proposal. 

As the final remark, we would like to comment that both the variables to evolve by the equations, and the variables 
used to confirm the Hermiticity of the system should be common throughout all evolutions. Otherwise, we cannot 
apply the energy inequality for the evolution of that system. From this point of view, we think it necessary to consider 
the secondary triad reality condition throughout evolution of this system. 

To summarize, we tried to follow ILR's procedure to remove the restriction of the triad reality condition in our 
system, which casts on our definition of symmetric hyperbolicity, and which is based on the fixed inner product as of 
its Hermitian form. We, however, see that ILR's procedure does not work in our system since it requires the restriction 
of the secondary reality conditions of the triad. Therefore we conclude that we cannot de-constrain restrictions any 
further. 
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[Additional Notes in this online version] 

Based on the results of this article, the authors developped a set of dynamical equations which forces the spacetime 
to evolve to the manifold that satisfies the constraint equations or the reality conditions or both as the attractor 
against perturbative errors. This report is available as: H. Shinkai and G. Yoneda, Phys. Rev. D60, 101502 (1999). 

The authors also performed numerical comparisons of three levels of the hyperbolic forms obtained in this article. 



The report is avilable as: H. Shinkai and G. Yoneda, ;r-qc/0005003 
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